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Abstract- Seeds of sunflowers are often modelled by the map n i — > 
<pe(n) = y/ne 2mnd leading to a roughly uniform repartition with two con- 
secutive seeds separated by the divergence angle 2tt9 for 9 the golden ra- 
tio. We associate to an arbitrary real divergence angle 2tt9 a geodesic path 
70 : M>o — > PSL2(Z)\H of the modular curve and use it for local descrip- 
tions of the image ye(N) of the phyllotactic map ipg. 

Given a real parameter 9, we call the map (pg : N — > C defined by 

the phyllotactic map of divergence angle 2ir6 (measured in radians). The 
image tpg(N) of a phyllotactic map is the phyllotactic set (of parameter 9 or 
divergence angle 2it9). A phyllotactic set tpg(N) is uniformly discrete (i.e. 
two distinct elements of ^e(N) are at distance at least e for some strictly 
positive e) with uniform density if 

9 = [a ; a%, a 2 , ■ ■ ■] = a + - — ^—^ — 



ai + 



0,2- 



is irrational with bounded coefficients oo, oi, 02, . • ■ in its continued fraction 
expansion. 

Among all possible parameters, the value given by the golden ratio 
= [1; 1, 1, 1, . . . ] (or closely related numbers) stands out and gives 
a particularly nice configuration. Figure 1 displays a few hundred small 
points of </? (1+v /5 )/2 (N). 

Finite approximations of ( / 9 (i+ v ^)/2(^) can be observed in capitula (heads) 
of sunflowers or daisies (the map fn+y/fyfti sometimes also called the sun- 
flower-map, has been proposed in [13] as a model for heads of sunflowers). 
Joining close points of tp^ / 2 (N) we get parastichy spirals appearing in 
pairs of crisscrossing families enumerated by two consecutive elements of the 
Fibonacci sequence 1, 2, 3, 5, 8, 13, 21, . . . . Explaining the occurence of the 



1 Keywords: Lattice, hyperbolic geometry, phyllotaxis, sunflower- map. Math, class: 
11H31, 52C15 Primary: 92B99 



1 



Figure 1: All points of (Pn+^n^) m ^ ne ^ sc i z e ^ I \ z \ — ^0}- 

golden ratio and of Fibonacci numbers in Botanies is the goal of Phyllotaxis, 
see for example Chapter XIV of [12] or [6j for more recent developments. 
The aim of this paper is to describe an elegant framework involving hy- 
perbolic geometry. The emergence of the golden ratio (or of closely 
related numbers) and of Fibonacci numbers enumerating families of paras- 
tichy spirals is then a consequence of natural constraints. 

Reasons for Phyllotaxis should be separated from the mechanisms which 
are involved. How Phyllotaxis works is surely best adressed by biologists, 
biochimists or biophysicists. The reason for Phyllotaxis is efficiency of some 
sort (a precise definition is perhaps not so easy) which can take several 
forms. It is perhaps a physical notion like energetic efficiency or it involves 
geometric quantities like isoperimetry (which leads probably ultimately also 
to some kind of energetic efficiency). The link between the two aspects is 
natural selection. To say it in a nutshell, ubiquity of phyllotaxis involves 
mathematics: A few geometric configurations optimize some natural quanti- 
ties. Thus they are favoured by living organisms through natural selection. 

Interestingly, the two aspects are spatially separated: Reasons for Phyl- 
lotaxy, due to a globally optimized quantity, are of an asymptotical nature. 
They are best adressed by studying the large part of a plant which is rela- 
tively far from the center consisting of the bud which is responsible for the 
growth-process of a flower. Asymptotic arguments are thus not a weakness 
but are relevant when trying to answer why phyllotaxis occurs (assuming 
Darwin's theory of evolution). 



An outline of the paper is as follows: We associate to a phyllotactic map 
(pg with real parameter 9 £ [0, 1) the curve 

Aiirt Voi^Qt 2 + 4wr(l - 9)t 
(0, +00) 9 * ^ 7 ,(t) = i -^ TT = 

of the Poincare halfplane i = {z G C ^(z) > 0}. Since ipg = t^e+n for all 
n E Z, we extend the definition of 79 to 6 E R by setting 

. , 4i7rt 
7«(*) 



4t7r{0}t + 1 

where 9 = [9\ + {9} with [9\ £ Z and {9} £ [0, 1) denoting the integral and 
fractional part of 9. It is easy to check that 70 defines a geodesic with respect 
to the usual hyperbolic metric of H. We call 70 the phyllotactic geodesic asso- 
ciated to ifg. For large N £ N, the projection of 70 (N) (often identified with 
70 (iV) in the sequel) onto the modular curve PSL2(Z)\IH classifying complex 
lattices up to similarity describes, up to an affine orientation-preserving sim- 
ilarity, the affine lattice obtained by "linearizing" the phyllotactic map ipg 
in a neighbourhood of tpg(N). Linearization is an asymptotical construction 
involving an error of order O (^y) m a neighbourhood of a point ip$(N). 
In particular, it breaks down for very small values of N. This failure should 
have no serious consequences: Indeed, the most interesting phyllotactic sets 
have obviously already good packing properties at their center. Moreover, 
the center yields a very small contribution to interesting quantities like en- 
ergy, mean isperimetric values of Voronoi domains etc.. 

The apparition of the golden ratio 1+ 2 V ^ can now be explained by the 
fact that the corresponding phyllotactic geodesic 1 1 — > rr^i ; avoids 

2iix\ — 1+v 5)t+l 

the cusp of the modular curve PSL2(Z)\H. This ensures small excentricity 
(or, equivalently, good isoperimetric constants) for Voronoi domains of the 
phyllotactic set 9?(i +v /5)/2(^) where the Voronoi domain V v of a point v £ S 
with respect to a discrete set S in a metric space E is the subset V v = 
{x £ E I d(x,v) = min^s d(x, y)} of all points closest to v. Plants having 
seed-areas with good isoperimetric constants should be favoured by natural 
selection since they need less material for constructing seed coats around 
seeds of given volume. Thus, Diophantine properties of the golden ratio 
r = ensure that ^(N) (or ipo(N) for 9 a close relative of the golden 

ratio t giving rise e.g. to the sequence 1, 3, 4, 7, 11, 18, . . . also observed in 
Phyllotaxis) satisfies natural constraints in the class of all sets of the form 

We can either consider that the points of the sequence ipo(0),ipg(l), . . . 
appear sequentially with constant divergence angle 2n9 between consecutive 
points on the so-called ontogenetic spiral 1 1 — > \/te 2lnet or, sticking perhaps 
closer to biological reality, we can consider a sequence of paths [— n, 0] 9 
t 1 — > s n fl(t) = \Jn + te 2l7Tne describing a situation where the n— th point 
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has appeared at the origin at time —n in the past. It has then slowly moved 
outwards on the halfray M>oe 2 " m6 ' until reaching its present location ipg(n) 
at the time t = 0. 

The factor y/n in the formula for <po(n) ensures that there are R 2 points 
(of roughly equal "importance" , measured for example by the area of Voronoi 
domains) in a disc of large radius R. Areas of Voronoi domains defined by 
</?e(N) are asymptotically equal to tt if 9 is irrational. 

Tools and part of the results of this paper can be adapted to a slightly 
more general situation given by functions t i — > p(t)e 2l7Tta ^ where p(t) is a 
suitable increasing function and where the local divergence angle 2-7ra(t) is 
allowed to vary very slowly. 

The sequel of the paper is organized as follows: 

Section [Q states the main result. 

Section [2] recalls a few well-known and useful facts concerning complex 
lattices and hyperbolic geometry. 

Section [3] contains identities involving continuous fraction expansions. 

Linearizations of phyllotactic sets are described in Section SJ 

We construct the phyllotactic geodesic 70 in Section This leads to a 
proof of Theorem 11.11 

Section describes a construction of a slightly different geodesic. 

Section [7] is devoted to metric properties of phyllotactic sets. 

Parastichy spirals are defined and studied in Section [8l 

Section [9] reviews briefly a few other models appearing in the literature. 

Finally, Section [10] discusses a possible experimental verification (or refu- 
tation) of the existence of a phyllotactic geodesic in real sunflower-capitula. 

1 Main result 

Let (E, dist) be a metric space, x an element of E and e, R two strictly 
positive real numbers. Two discrete subsets A, B of E are e— close in the 
open ball of radius R and center x if there exists a map tp : A' — > B' which 
is one-to-one and onto between subsets A' C A and B' C B containing all 
points in A, respectively B, at distance at most R from x and which moves 
all points of A' by less than e, i.e., we have dist(a, VK a )) < e for all a £ A'. 
Intuitively, two discrete sets A, B are e— close in the open ball of radius R 
centered at x if A and B are "equal up to an error of e" in (a neighbourhood 
of) the ball of radius R centered at x. 

Theorem 1.1. Given e > and R > 0, there exists an integer N = N(e, R) 
such that for every 9 £ [0, 1) and for every n > N , the set (pg(N) is e— close 
in the open disc of radius R centered at ipo(n) to an affine lattice in the 
equivalence class (i.e. up to orientation-preserving affine similarities) of 
Z + Z . 
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The map 



(0,oo) 3 t 



le{t) 



1 + 4m{9}t 



4lTTt 



16n 2 {9}t 2 + Amt 

1 + {^{9}t) 2 



defines a geodesic of the hyperbolic Poincare halfplane (endowed with the 
hyperbolic metric — at z = x + iy G H) , see for example Lemma 12.11 We 



Remark 1.2. (i) Theorem li.il holds for rational 9: In this case the phyl- 
lotactic geodesic ends up in the cusp of the modular curve PSL2(Z)\HL The 
associated affine lattices A# n degenerate into discrete subgroups of rank 1 
in the sense that they intersect a ball of fixed radius R centered at an affine 
lattice point along a translated copy of a discrete subgroup having rank 1. 

(ii) Denoting by V(n) the Voronoi domain of (fo(n) G </_>(N) we have 
lim n _> 00 volV(n) = 7r if 9 is irrational. 

For rational 6 = | with with q > 3 and p, q coprime integers, the union 
of Voronoi domains (defined by <pg(N)) of all points at distance < R from 
the origin is essentially a regular polygon with q sides and inradius R. This 

implies lim.n_j.oo volV(n) = gtan ^ = 7r + J^p + O (^^j ■ 

For 9 G all Voronoi domains are unbounded and thus of infinite 

volume. 

The existence of phyllotactic geodesies gives a measure of "similarity" 
of phyllotactic sets in neighbourhoods of <pe{n) and (pe'{m) by considering 
the hyperbolic distance _m(7„ PSL/2(Z)7e/(m)) between the two orbits 
PSL2(Z)7e(?7,) and PSL2(Z)7g/(m). Statement (ii) of Remark 11.21 ensures 
that a small distance de(7e(n), PSL2(Z)7g/(m)) implies the existence of bi- 
jections between and tpg'(N) which are almost isometries in neigh- 

bourhoods of (fo(n) and <pQi(m) if 9 and 9' are irrational and have continued 
fraction-expansions with bounded coefficients. More generally, this holds if 
70 (n) (and thus also 70' (m)) is far from the cusp of PSL2(Z)\BL 

2 Complex lattices and hyperbolic geometry 

For the convenience of the reader, we recall a few elementary and well-known 
facts first of the theory of lattices, following closely parts of Section 2.2 in 
Chapter VII of [11 j . then of hyperbolic geometry, see for example [1]. 

2.1 Lattices of C 

A lattice in C is a free additive subgroup generated by two M— linearly in- 
dependent elements wi,W2 of C. In the sequel, we consider lattices only up 
to orientation-preserving similarities. Two lattices F and A of C are thus 




70 (t) has 
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equivalent if A = Ar for some non-zero constant A G C*. Given a basis 
LJi, 0J2 of a lattice T = 7Luj\ + Zto>2, we consider z = Up to replacing, say, 
u>i by — cji, we can suppose that the imaginary part y = Q(z) of z = x + iy 
is strictly positive. 

Thus, a lattice Zu>i + Z^2 — fz + ^z^j ^2 is equivalent to the lattice 
T(z) = 'L + 'Lz generated by 1 and by the element z = — of the open upper 
half-plane H = {z G C | S(z) > 0}. 

Given an unimodular integral matrix 5 = f ^ J G SL2(Z), the quo- 
tient 



z 



u>2 

oj'r. cuj\ + duo c— + d 

UJ2 



associated to the basis uj[ = aoj\ + 6^2 , ^4 = cw i + ^2 of a lattice Zu>i + Zw2 
is obtained from z = ^ by the usual action 5.2; = of the modular group 
PSL 2 (Z) = SL 2 (Z)/ ± Id on H. 

Hence the map "Luii + Zil>2 i — > z = ^ G EI induces a one-to-one cor- 
respondence between equivalence classes <C*r of lattices and points of the 
modular curve PSL 2 (Z)\H, see Chapter VII, Proposition 2 and Proposi- 
tion 3 of [llj. A fundamental domain for the action of the modular group 
PSL2(Z) on H is given by the fundamental domain 

M = iz G C I \z\ > 1 and |»(z)| < ij (1) 

for PSL2(Z)\H. Two elements 21,22 of Ai represent the same equivalence- 
class of lattices if and only if either Z2 = — — or z^ = Z\ ± 1 . The modular 
curve PSL2(Z)\H = A4/ ~ is a complex orbifold with two conical points 
represented by -1 + tv/ ^ ( f angle ^ and corresponding to regular hexagonal 
lattices) and by i (of angle tt and corresponding to square lattices) and 
with a cusp (corresponding to a neighbourhood of the degenerate case of an 
additive subgroup of rank 1 in C). 

An affine lattice is a coset a + T obtained by translating a complex 
lattice r C C by some vector a G C. We consider affine lattices only 
up to orientation-preserving affine similarities. Equivalence classes of affine 
lattices are also in one-to-one correspondence with elements of the modular 
curve PSL 2 (Z)\M. 



2.2 Hyperbolic geometry on the Poincare half-plane 

We recall a few facts concerning the hyperbolic Poincare half-plane EI, see 
[T] for an elementary introduction to hyperbolic geometry. 

The upper half-plane EI = {z G C \$s(z) > 0} can be turned into a real 
hyperbolic simply connected Riemannian manifold of dimension 2 and of 
constant curvature —1 by equipping it with the Riemannian metric (ds) 2 = 
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dx at a point z = x + iy £ M. The Poincare half-plane is the hyperbolic 
manifold (still denoted by) H obtained in this way. 

The group of all orientation-preserving isometries of the Poincare half- 
plane is given by the set of all Mobius transformations 

a b \ az + b 
c d J cz + d 

defined by matrices ( ^ j in SL2(R), respresenting elements in PSL2(R) = 
SL 2 (K)/±Id. 

The geodesies of EI are half-circles (with respect to the usual Euclidean 
metric of C) centered at the boundary R of EI C C or halnines {a + iy € C | 
y > 0} C EI perpendicular to R. 

An orientation-preserving isometry t of the Poincare half-plane is hyper- 
bolic if it admits an invariant geodesic on which it acts by a translation. A 

a b 



Mobius transformation associated to I d J ^ ^L 2 (R) defines a hyper- 
bolic isometry l if and only if |a + d| > 2. The invariant geodesic of i is given 
by the halfcircle in EI delimited by the two real points d ±V(^ a) +4ftc 
c ^ respectively by the halnine {b/(d — a) + iy \ y > 0} otherwise. 

Lemma 2.1. For f ° ^ j G GL 2 (R) with cd ^ and positive determinant 
ad — be > 0, i/ie image of the map from R>o = (0, +oo) into C defined by 

^ ait + b 

cit + d 

is an open half-circle ofW centered on a< ^ c with radius I a< ^ c I (u>i£/i respect 
to the Euclidean norm \\ z \\= \Jx 2 + y 2 for z = x + iy G C). 

Proof Mobius transformations preserve geodesies of H. Thus the Mobius 

a b 



transformation defined by the matrix I ^ J sends the geodesic {ti \ t > 

0} C EI onto a geodesic of EI with finite boundary points given by ^ (corre- 
sponding to t = 0) and | (corresponding to t = +oo). This geodesic is the 
open halfcircle (for the usual Euclidean metric of C) of the upper halfplane 
with center \ (§ + f) = and diameter |f - a^s| = |^^|. □ 

Remark 2.2. Lemma \2.1\ is equivalent to the identity 

( act 2 + bd ad + bc\ 2 . / ad — be \ 2 2 f ad — be 



+ 



c 2 t 2 + d 2 2cd / \c 2 i 2 + d 2 / I 2cd 
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3 Continued fractions 



We denote by [ao; a±, 02, ■ ■ ■ ] the continued fraction expansion 

= a H — ; — , 

ai H i — 

a2+ ' • 

of a real number 0. The coefficients ao,a\,... are recursively defined by 

<H = \ 0i \ where On = and n = „ — = ffl 1 , if n -\ 4 Z, respectively 

by #n = if # n _i G Z. The coefficient ao of a continued fraction expansion 
can be an arbitrary integer (positive, zero or negative). 01,02, .. . are either 
all strictly positive or they start with a finite number of strictly positive in- 
tegers followed by an infinite string of zeros. The last case arises if and only 
if is rational. The sequence ao,oi, • • • is infinite and ultimately periodic 
with non-zero period if and only if Q[0] is a quadratic number field. Every 
irrational number has a unique continued fraction expansion. Rational num- 
bers have two expansions given by [ao; a±, . . . , o m , 1] and [ao; 01, . . . , a m + 1] 
for suitable integers m > 0, oq G Z, a\, . . . , a m > 1. 
We have the continued fraction expansions 

n = [a n \ a n +i, a n+ 2, a n+ 3, • • • ] (2) 
for all n 6 N. Convergents for 



io = a -\ r 

ai H - 



Ct2- 

are rational numbers of the form 

P-2 _ p_l _ 1 fc_ + QnPn-1 

9-2 1' 0'on g n -2 + a n q n -i 



[a ;ai,...,a n ], n > 



and can also be define by 2a = [a - Ql . . . a 1 see Theorem 149 in HI. 
Intermediate convergents are given by 

Pn-2 + fcgn-l , _ /n 

— , fc 6 {0, . . . , a n - 1} . 

On_ 2 + kq n -i 

The easy identity 

Pn-lQn ~ PnQn-1 = (-1)" (3) 

(see Theorem 150 in HJ or Theorem 2 in [7]), equivalent to ^— - — — = 

(-1)" 



<?n-l<?i 



-, implies 



P2n ^ P2n+2 ^ ^ a ^ ^ P2n+1 . P2n-1 
< < • • • < a < • • • < < 

l2n l2n+2 <l2n+l Q2n-1 
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This shows 



Pn 
Qn 



< 



(4) 



and ensures that convergents are excellent rational approximations of an 
irrational number. 

The following result is essentially identity 10.3.2 of [I]: 



Proposition 3.1. We have 



for all n > 0. 

Proof The result holds for n 
We have 

Pn-1 + On+lPn 



Pn-2 + OnPn-1 
Qn-2 + QnQn-l 



0. 



Pn-1 + 



(Pn-2 + ttnPn-l' 



-1 + On+lQn 



Qn-1 + e„-a n (Qn-2 + a n q n -i) 
Pn-2 + OnPn-l 



Qn-2 + n q n -\ 
which ends the proof by induction. 
Lemma 3.2. We have 



□ 



Qn-1 + Qn+lQn 



ln-2 + v n q n -l 



Proof We have 



>n+l 



Qn-1 + Vn+iq n 



ir^n-l + Qn 
1 

9?i — a n )g n _i + q n 
1 



VnQn-l — anQn-1 + 9n-2 + fln^n-l 

where we have used the recursive definitions = 9 n — a n and g n 

g- n _ 2 + fln^n-i of 6 n+1 and of g n . 

Proposition 3.3. We have 

p n -2 + Xp n -i 



□ 



"nQn—l) 



(5) 
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Proof Proposition 13,11 shows that the result holds for x = 6 n . 
Since ([5]) is equivalent to the identity 

9 n — x 

(g n _ 2 + Xq n -i)9 - (p n -2 + Xp n -i) = V — (-1)™ (6) 

Qn-2 + fnQn-1 

involving affine functions of x, it is enough to show the equality 

Qn-10 ~ Pn-l = — \> • (7) 

Qn-2 + "nQn-l 

This holds for n = since it boils down to — 1 = —1. By induction, we have 
for x = a n the identity 

Pn-2+dnPn-! (6 n - a n )(-l) n 



(8) 



Qn-2 + &nQn-\ (Qn-2 + UnQn-l) (Qn-2 + OnQn-l) 

which can be rewritten as 

e Pn_ (-IT 

Qn Qn+lQn(Qn~2 + OnQn-l) 

using the recursive definitions of p n , q n and n +i. 
The identity 

&n+l (Qn-2 + OnQn-l) = Qn-l + &n+lQn 

equivalent to Lemma 13.21 yields now ([7|) for n + 1. □ 

Remark 3.4. Identity (corresponding to the specialization x = a n of 
Proposition 1 3. 3j) strengthens inequality ^ since a n+ \ = \ 6 n+ \\ < 9 n+ \ and 

Qn = Qn-2 + O-nQn-X < <7n-2 + &nQn-l- 

Lemma 3.5. We /iawe 

1 On+l 1 



+ "nQn—1 Qn—l + dn+lQn Qn + #n.+29n+l 

Proof Using the identities 

Qk = Qk-2 + a>kQk-i, 
Ok+i = - a , 

"k — a k 
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for k = n and k = n + 1 we have 

1 CLn+l 



+ "nQn—l Qn— 1 + Qn+lQn 

1 a n+ i 



In — o, n q n -\ + + Gn+llr, 

1 a-n+1 

+ JTT^n-l Qn-1 + On+lQn 

u n-\-l 

#n+l ~~ a n+l 
Qn-1 + 9n+lQn 

Gn+l — Qn+1 

+ fln+l^n + {0 n +l ~ 0"n+l)Qn 
1 



#n+2<7n+l + Qn 

which ends the proof. □ 

4 Linearization 

Proposition 4.1. We /iaue 

+%(o,6) 

where q^ is the denominator of the k—th convergent 2fc = [oo; ai, ■ ■ ■ , afc] /or 
= [ao; ai, 02, . . . ], where Ok = [a^; ajfc+i, a/c+2> • • • ] *s defined as in HJ) and 
where the error Ej(a,b) is asymptotically given by 

(i^f (aqj-i + bq,) 2 + ^=K-l + b Qj )5 - 2vr 2 ^ 2 ) e M 
6= -a ( - 1)J >' +6- (-'^ 



if \ip 6 (n + aq j _ 1 + bqj) - y?e(n)| = 0(1). 

Corollary 4.2. If qj-i and qj are denominators of two consecutive conver- 
and^ 



gents |— !- and ^ of 9 such that qj-i < \/n < qj, then the smallest points 



of 
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are close to the smallest points of the lattice 
z ( «t± _ { _ iy Wn» i+1 ) +Z (^ + { _ iy 2^i + 2 \ (g) 

with an error of order O (^~~^J ■ 

The lattice described by Q contains always a non-zero element of abso- 
lute value smaller than \J\ + 47r 2 < 3n. 

Remark 4.3. Fundamental domains of the lattice A given by ^) have area 
it as shown by the identities 



(-l) J 'det 



2y/n qj-i+dj+iqj 

qj (-l)i2iwy/nO j+ 2 
2\/n <lj+8i+2<lj+i 



Vj+2 Vj+l 
^TTTfl a + Qj a TO 7T n 

1 ^i+i 



7T 



where the second equality is given by Lemma \3.?i 

Since the regular hexagonal lattice has maximal density, the lattice A 
contains always a non-zero element of absolute value at most ~ 1.9046. 

Proof of Proposition 14.11 Setting 

F(s, 7) = Vn~+~s e 2i7Tnd+2i ^ 

we want to approximate 

F(aq j . 1 +bq j ,5)-F(0,0) 

where 5 = (aqj-i + bqj)9 — c is the difference between (aqj-i + bqj)6 and the 
integer c closest to (aqj-i+bqj)9. We do this in the usual way by considering 
the linear approximation 

OF dF 
L = — (0,0)(a<&_i + bqj) + ^-(0,0)5 

and by estimating the error using second-order derivatives. The necessary 
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partial derivatives of F are: 

dF 1 



d 2 F 



dF 

(0, 0) = 2ivrVne- 



2iirn6 



2inn9 



— (0 0) = ~ l e 2i ™ 6 



m 2mn6 



(0,0) = — e l 
(0,0) = -47T 2 v / ^e 2 ^ n9 . 



dsd^f ' \/n 
d 2 F 
c?7 2 



The contribution coming from ^-(0, 0)(aqj-i + bqj) to L is given by 

1 



In order to compute ^(0,0)5 we split 5 into 5 = a8j—\ + bSj where S/. for 
k £ {j — is the difference between q^O and the integer closest to q^d. 
Since % is a denominator of the convergent = [ ao ; ai, . . . , a^l of 0, this 
integer is given by the numerator p/,. We have 

a Pk Pk + 0k+2Pk+l Pk 



Qk Qk + #fc+2<7fc+l qk 

(Pk+lQk — PkQk+l)^k+2 



Qk(.Qk + 8k+2Qk+l 



dk+2 



qk(qk + 6k+2Qk+i) 
where we have used Proposition 13.11 and identity ([3]). This yields 

9k+2 



(-l) k 



Qk + ^k+2Qk+l 



and shows 

dF dF 



(0,0)5 = — (0,0)(a<5 i _i + bSj) 



ds V Qj-i + Oj+iQj Qj + 8j+2qj+iJ 
The order of the error is given by 

5 1^(0. + b Qj ) 2 + J^(0, 0)(ag J _ 1 + bqj)5 + ±0(0, 0)5 2 

and can be evaluated easily. □ 
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5 The phyllotactic geodesic and proof of Theorem 

Using (|3.2p we can rewrite the lattice A given by formula ([9]) of Corollary 



as 

201 qj-2 + Ojqj-iJ \2^/n qj-i + 6 j+1 qj 

In particular, the lattice A is similar to the lattice Z + Zrj(n) where 

._, _ ii7r| 

-iW = -(-^ 7^1 ^ ■ (10) 

Remark 5.1. A straightforward computation shows that the imaginary part 
ofTj(t), given by 

<?. + igggg a ' 

is strictly positive if t is strictly positive. 
Theorem 5.2. For all j > we have 

(r j (t)-(-iya j+1 )T j+1 (t) = -l . 
Proof Theorem 15.21 boils down to the identity 

, .• Aiirt • 4iirt 

qj-2 + Vjqj-i qj-i + Vj+iqj 

■ Aiixt 

= q 3+ i-{-iy—— — - — . 

Qj + Vj+2qj+l 

The identity qj+i = qj-\+aj + \qj shows that the constant parts (with respect 
to t) of both sides are equal. Linear coefficients of t are equal by Lemma 
1331 □ 
Proof of Theorem 11.11 By Theorem 15.21 the two geodesies defined by r,- 
and Tj+i are related by the integral Mobius transformations 

Tj+1 = Tj - (-i)'Vx = ( -°i (-iy aj+1 ) • T > 

and 



(-iy aj+1 T j+1 - i = i' (-iy aj+] 

Tj+l 



■T j+l . (12) 
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Thus they project onto a unique geodesic on the modular curve PSL2(Z)\H 
represented for example by 

(where = 1, 5-1 = 0, g = 1 and 9\ = = jgj). 

This implies Theorem 1 1 , 1 1 since the linearization error is of order O {^^j 
for elements of <po(N) at bounded distance from <pg(n). □ 

Formula ()13p defines a geodesic of the hyperbolic half-plane for every real 
number 9. Indeed, f)13|) is a vertical half-line (and thus a geodesic) if is 
integral and it defines a halfcircle of H orthogonal to R (and thus a geodesic) 
with boundary points corresponding to t = and corresponding to 
t = oo otherwise. 

Diophantine properties of 9 are related to the dynamical behaviour of 
the geodesic 70 projected onto PSL2(Z)\H as follows: after starting at the 
cusp, (the projection of) 70 turns (slightly less) than a\ times around the 
cusp before passing between the two conical points of the modular curve. 
It turns then in the same sense (and slightly less than) ai times around 
the cusp before crossing again the shortest geodesic segment joining the two 
conical points and so on. A large coefficient au causes the (projection of the) 
geodesic 70 to climb the modular curve up to a height given asymptotically 
(in ak) by a^i G M. This gives rise to points of ipe(N) having Voronoi 
domains with bad isoperimetric properties. 

For a divergence angle 2ir9 determined by the golden ratio 9 = (or 
close relatives of it) the continued fraction expansion involves only ones (or 
only ones after perhaps a few initial "accidents"). This is the optimal situ- 
ation leaving no possibility of improvement. In particular, the phyllotactic 
geodesic 7( 1+v ^)/ 2 * s asymptotically equal to the geodesic 

(-1 + V5)it - 1 - y/h _ -1 - y/E+ (-1 + V5)t 2 + 2ity/5 
1 1— »• 7W - 2(it + l) " 2(1 + i 2 ) 

with boundary points ~ 1 ^ v/ ^ (and containing the points = 7 ( 

and % = 7 O^^j of M). The equality 

l(t) + 1 



7(*) + 2 



shows that 7 is invariant under the integral Mobius transformation defined 
by the matrix f j 2 ) = ( 1 1 J ' ^ projects onto the shortest closed 
geodesic of the modular curve PSL2(Z)\H. 
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Since ^j{t) = 7^w> the parametrized geodesic 7(f) has the same in- 
stant speed QfL t \\ = j as the phyllotactic geodesic 70 (i). We have 
thus asymptotically 79 (<) ~ 70 (0 t) . The phyllotactic set ipg (N) has thus 
almost isometrical neighbourhoods around <pg(n) and (pg(m) if n is large and 
m is close to 9 A n. 



6 A geometric construction 

We construct in this section a slightly different geodesic on the modular 
domain PSL2(Z)\KI which is asymptotically associated to linearizations of 
W(N). 

We denote by £ = Cg the line R(l, —9) of slope —9 containing the 
origin. A convergent ^7 = [ao; cti, ■ . . , a,j] of 9 = [ao; ai, • • • ] yields an 
integral point (qj,—pj) close to £ as follows: Denoting by ftc{qj-,—Pj) the 
orthogonal projection of (g^, — Pj) onto £ we have 

-P,) = -0) + O ^ (1, -0) 

for all j > 0. The equality qj9 — pj = (~l) J ' g . + g J ^ 2 2 g . +1 corresponding to the 
case n = j + 2 and x = of Proposition 13.31 implies the identity 



(<&, -Pi) -<&(!, -0) = 0,(-l) 



The orthogonal projection Tt c ±(qj,—pj) of (qj,—pj) onto the line £ J " = 
R(#, 1) orthogonal to £ = R(l, — 9) is given by 

(-iyg 3 - + 2 (M) 

(<& + 0i +2 g j+ i) (1 + 02) • 
Thus we can rewrite the right side of the obvious identity 

(lj,-Pj) = Tc(qj,-Pj) +n c ±{qj,-pj) 

as 



We endow now Mr with an Euclidean metric dsg n turning the vectors 
(e.i) . 

' 47rn.(l+0 2 ) 



(1, —6), 4 (ilgi\ into an orthogonal basis. Comparision of (fTi]) with ([9]) 
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shows that the lattice (Z 2 ,dsg jn ) is asymptotically equivalent with the lin- 
earization of ipo (N) at the point (po(n). Since 

(M) = -^(l,-e) + 47rn6 



(9 A) 



(M) = -TT^2^-°)+^ n 



l + 9 2K ' ; 47rn(l + 2 )' 

the Euclidean lattice (Z 2 ,dse jn ) corresponds to the point of PSL2(Z)\IH rep- 
resented by 

1 + <Hnn9(l + B 2 ) 
9 - Aiimil + 2 ) 

9 - 16ir 2 n 2 6(l + 6 2 ) 2 + 4i7rn(l + 6 2 ) 2 
9 2 + 16Tr 2 n 2 {l + 9 2 ) 2 ' 

All these points are elements of the hyperbolic geodesic with boundary points 
| (for n = 0) and — 9 = ( °° q ^ {IJ ^ or n = °°) - 

Remark 6.1. Since we have 

1 + 4iTTt9(l + 9 2 ) ( a 1 \ 4i7rf 



6> - 4z7rt(l + 2 ) V- 1 / 4w{^}t + l 
1 + 4i7rf6»(l + 6» 2 ) Amdt + 1 



+ 



6> - 4ivrt(l + 6> 2 ) 4ivrt 

B_ 

Amt{9 - Amt{\ + 9 2 )) 

1 + 9 2 9 2 



9 2 + 16vr 2 t 2 (l + 9 2 ) 2 4vrt {9 2 + 16vr 2 f 2 (l + 2 ) 2 ) 
the hyperbolic distance between the two points 

1 + 4i7rtfl(l + 2 ) / a 1 \ 4i7rt 
- 4i7rt(l + 2 ) V — 1 / 4i7r{fi}t + l 

o/H is roughly given by 



? 2 + 16vr 2 t 2 (l + 2 ) 2 



1 + 2 2 



47rt(l + 2 ) 2 
which simplifies to 

1 fl2 



2 + 16vr 2 i 2 (l + # 2 ) 2 47rt (6> 2 + 16vr 2 i 2 (l + 9 2 ) 2 ) 



4?rt(l + # 2 ) 167r 2 t 2 
T/ius ii is asymptotically equal to 47rf ^ + g^) which is asymptotically much 
smaller than the error O (~^=\ due to linearization at a point ipe(n) = 
^h~e 2i7Ten of order 0{y/i). 
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7 Metric properties of phyllotactic sets 



A subset S of a metric space E is uniformly discrete if there exists a strictly 
positive real constant 5 such that d(a, b) > 5 for every pair a, b of distinct 
points in S. Equivalently, S is uniformly discrete if open balls of radius 
5/2 centered at all elements of E are disjoint (for a small strictly positive 
constant 5). 

A subset S of a metric space E is an e—net if every point of E is at 
distance at most e from a point of S. Equivalently, E is covered by the set 
of closed balls of radius e centered at elements of S. 

The following result is a straightforward consequence of the fact that 
areae of Voronoi domains defined by ipg (N) are asymptotically equal to it if 
6 is irrational: 

Proposition 7.1. The following assertions are equivalent: 

(i) 6 is irrational and has bounded coefficients a±, a,2, ■ ■ ■ in its continued 
fraction expansion 6 = [do; a,\, a,2, ■ ■ ■ ]■ 

(ii) </?£)(N) is uniformly discrete in C (identified with the Euclidean plane 
in the obvious way). 

(Hi) <po(^) is an e—net of C. 

(iv) All Voronoi domains of tpg(N) have bounded diameter. 

(v) Discs of radius R (and arbitrary centers) in C contain R(R + O(l)) 
points of ipg (N) . 

(vi) The image Jg([l, oo)) of the phyllotactic "half- geodesic" is contained 
in a compact subset of the modular curve PSL2(Z)\H. 

We leave the proof to the reader. □ 
Observe that ipg (N) is never uniformly discrete if 9 is rational. 

8 Parastichy spirals 

We denote by dA4 the boundary in C of the fundamental domain M defined 
by ([!])• The interior M \ dM corresponds to lattices having a unique pair 
±n of opposite shortest non-zero vectors and a unique pair ±v of shortest 
vectors which are M— linearly independent from ±u. More precisely, for 
z £ A4 \ dA4 the unique pair ±u of non-zero shortest vectors in Z + Zz is 
given by ±1 and the unique pair ±v of shortest vectors outside M coincides 
with ±z. Notice that IR— linear independency of v from u is necessary in 
order to discard ±2n, ±3n, . . . which might be smaller than v for lattices 
associated to z G M. with large modulus. 

Lattices corresponding to elements z of norm 1 in Ai have (at least) 
two pairs of shortest vectors given by ±1 and ±z in Z + Zz. The regular 
hexagonal lattice corresponding to the unique lattice with three 

pairs ± 1, ± 1+ ^ v ^ , jr 1 ~2 °f shortest non-zero vectors. Lattices associated 
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to z = ^y^- for t > y/3 have a unique pair ±1 of shortest non-zero- vectors 
and two non-real pairs ±z and ±(z — 1) of shortest non-real vectors. 

Connecting points of a lattice A indexed by z S M. \dM. with their closest 
neighbours we get a set of parallel lines. Joining closest lattice-points on two 
such adjacent lines we obtain a second set of parallel lines. These two sets of 
parallel lines cut the complex plane into fundamental domains for A given 
by isometric rhombi. 

Parastichy spirals are analogues of these lines in ipg(N). More pre- 
cisely, we define (generically) the primary parastichy spirals of (fo (N) as the 
piecewise-wise linear paths obtained by joining vertices of ^(N) to their 
two approximatively opposite nearest neighbours. Similarly, we construct 
secondary parastichy spirals by joining vertices of <p$(N) to their nearest 
neighbours on adjacent neighbouring primary parastichy spirals. 

Primary parastichy spirals exist essentially at every point far from the 
origin except where they become blurred with secondary parastichy spirals. 
At such points (corresponding to crossings of the phyllotactic geodesic with 
the image of the unit circle in PSL2(Z)\H), primary and secondary paras- 
tichy spirals get exchanged. We call such a situation a parastichy transition 
of type I. 

Secondary parastichy spirals are however well-defined only if the local 
situation corresponds to a lattice indexed by an element of A4 which is not 
too close to the cusp. For example, the phyllotactic set p$(N) associated to a 
rational number 8 = | is contained in q half-rays originating at 0. Far from 
the origin, primary parastichies (and are no longer spirals) coincide with 
these half-rays and secondary parastichies make no longer sense. Moreover, 
for points z £ Ai with real part close to 1/2 (or —1/2) a family of secondary 
parastichy spirals fades away and is replaced by a new family of secondary 
parastichy spirals, giving rise to a parastichy transition of type II. A coeffi- 
cient en > 1 yields aj — 1 parastichy transitions of type II. The occurence of 
parastichy transitions of type II is easy to detect visually: it leads to much 
less uniform point distributions in (pg(N). Figure 4 displays two examples. 

The geometric construction of Section [6] shows the well-known fact that 
primary and secondary parastichies form two sets of spirals with different 
orientations if 9 is irrational. Indeed, primary, respectively secondary, paras- 
tichies around tpe(n) are defined by (pg(n±qj), respectively by tpg(n±(qj—i+ 
kqj)) for suitable integers j, k, see Figure 2 where one has to think of (0,0) 
as the point <po(n) and of C as the ray defined by M>o<£>0(n). (As always, 
qj—l and qj are denominators of convergents for 6.) The same integers j,k 
work for all n in some interval of large length compared to y/n. The plane 
M 2 is increasingly squezed (for increasing n) in the direction of C and ex- 
panded in the orthogonal direction C in the construction of Section [6l This 
implies that parastichies of both kinds bend away from the rays issued by 
the origin. Moreover, there is exactly one parastichy family of larger, re- 
spectively of smaller slope than <pg(n) as can be seen by inspecting Figure 
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2. This explains the apparition of crisscrossing spirals in Figure 1. Sec- 
ondary parastichy families are however no longer discernible (to the eye) 
for larger values of k, see Figure 4 where there are regions without obvious 
secondary parastichies. For irrational 9, they can however always be drawn 
"by continuity" , if we start in suitable regions where no problems occur and 
if we push them forward using type II transition for (the projection of) 
crossing the infinite boundary segment of Ai. 

8.1 Transitions for parastichy families 

The number of "parallel" primary parastichy spirals forming a common fam- 
ily is always a denominator qj of a convergent ^ = [°o> °i> • • • > a j] of the di- 
vergence angle 2tt0 = 27r[ao, a±, a,2, ■ ■ ■ ]■ The number of secondary parastichy 
spirals in a common family is a denominator qj + fcffj+i, k £ {0, 1, . . . , aj+2} 
of an intermediate convergent. 

The exterior region of Figure 1 for example contains 55 primary paras- 
tichy spirals turning clockwise and 34 secondary parastichy spirals turning 
count er clo ckwise . 

The evolution of the numbers of parastichy spirals (PS in the following 
table) can be described by: 











/ 


prim. PS 

qj-i sec. PS 


|7o(n)| ~ 1 


qj sec. PS 
qj prim. PS 


II 


qj-i + qj sec. PS 


|»(7fl(")l ~ 1 


qj prim. PS 


II 


qj-i + 2qj sec. PS 


m-Yo(n)\ ~ i 


qj prim. PS 


II 








/ 


qj-i + a j+1 qj = q j+1 sec. PS 
qj + \ prim. PS 


\le{n)\ ~ 1 


qj prim. PS 
qj sec. PS 


II 


qj + i prim. PS 


mie(n)\ ~ i 


qj + qj + i sec. PS 











where jg (n) denotes a representant of of 70 (n) in A4. 

Figure 2 attempts to illustrate the occurence of parastichy transitions 
using the geometric construction. 

Figure 3 shows the death and birth of a family of secondary parastichy 
spirals corresponding to a value of n such that |K(7#(n)| ~ ^. The verti- 
cal arrow represents a ray issued from the origin. Primary parastichies are 
represented by fat segments, the dying family of secondary parastichies is 
drawn with dotted segments and the newborn family of secondary paras- 
tichies is given by ordinary segments. Turning around the origin on a circle 
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(qj+i,-Pj+i) = A + a j+1 B 




Figure 2: Parastichy transitions 




Figure 3: Death and birth of a family of secondary parastichy spirals 

of constant radius ~ y/n and counting the numbers a of primary parastichy 
spirals, b of dying secondary parastichies and c of newborn secondary paras- 
tichies, one gets the relation c = a + b. More precisely, the piecewise linear 
path involving only segments of primary and secondary dying parastichies 
giving the best approximation of the circle with radius y/n consists of b seg- 
ments on primary parastichies and a segments on secondary parastichies. In 
order to work with the family of newborn parastichies, one has to replace 
every segment of the dying family by two segments, one from a primary 
parastichy and one from a newborn secondary parastichy. The number of 
segments on secondary parastichies (which is equal to the number of curves 
in the primary family) remains thus constant (and equals a) and the number 
c of segments on primary parastichies (which equals the number c of curves 
in the newborn secondary parastichy family) increases by a to c = a + b. 

In order to prove that parastichy families are enumerated by denomi- 
nators of (intermediate) convergents, it is now enough to remark that the 
assertion holds for the final number of primary parastichies if 6 is rational. 
A continuity argument implies the result in general. 
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Figure 4: Small points of (po(N) for 
(right side). 



1765- 



2858 



^ (left side) and 



The left half of Figure 4 shows all points of (p$(N) with real and imaginary 
parts smaller than 30 for = l7f %^ and 9 = e -1 . We have ^^"/^ = 
[0; 1, 1, 1, 1, 1, 1, 3, 1,1,1,1,...] and the first convergents are 

1 1 2 3 5 8 29 37 



1' 1' 2' 3' 5' 8' 13' 47' 60'"' 

The 13 parastichy spirals corresponding to the denominator 13 are clearly 
visible. This is of course due to the "large" coefficient 3 in the continued 
fraction expansion of 9 which leads to three parastichy transitions of type 
II. 

Similarly, we get for the parameter 9 = e~ l = [0; 2, 1, 2, 1, 1,4, 1, 1,6, 1, 1, 8] 
corresponding to the right half of Figure 4 the convergents 

1 1 3 4 7 32 39 



1' 2' 3' 8' 11' 19' 87' 106' 



with a clearly visible parastichy family corresponding to the denominator 
19. This family is due to the large coefficient 4 (the family corresponding to 
6 becomes visible at a larger scale) which causes four parastichy transitions 
of type II. 



8.2 Monodromy 

One can consider two notions of monodromy for phyllotactic sets: 

A first notion consists in moving a chosen basis for the local "lattice" by 

comparing bases of close points in the obvious way. No monodromy arises 

in this way: A closed loop gives rise to the identity. 

A second, slightly more interesting feature is translational monodromy: 

going counterclockwise around the origin on a piecewise linear path and 

stitching the obtained lattice elements (with respect to "bases" which are 
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"moved" continuously) together, we get at a point a + ib G ^(N) a vector 
close to 2ir(—b + ia). More precisely, using a basis V\,Vi associated to 
primary and secondary parastichy spirals, this vector is of the form (qj—i + 
kqj)V\ ± q.jV2 (with signs depending on the sign conventions for V\ and V2) 
if (fo(N) contains qj primary and (qj-i + kqj) secondary parastichy families 
at distance V« 2 + b 2 from the origin. 



9 Other models 

9.1 Cylindric and logarithmic models 

In [3] Coxeter, following [2], models the structure of pineapples or pine- 
cones by approximating their shape with a cylinder which he develops on 
the plane thus getting an infinite strip of a lattice. The visible features 
(scales) on pineapples are the Voronoi domains of this lattice. The obtained 
lattice should stay close to the hexagonal lattice which has optimal packing 
and covering properties. Working with a cylinder of circumference 2ir, we 
have thus to choose the optimal divergence angle 2tt9 such that the complex 
lattice 7L + Z(0 + ei) is close to the hexagonal lattice for small e. There is 
no exact control over the value of e since pineapples or pine cones are not 
exact cylinders. We should thus choose the value of such that the geodesic 
1 1 — > 9 + ti is overall optimal for small positive e. For a fixed value of 9, the 
map e 1 — > 9 + ei defines again a geodesic of the hyperbolic half-plane. 

The best choice is of course again given by A = 1+ ^ 2tt (or by its neg- 
ative) modulo 2tt yielding a geodesic which is asymptotically close to the 
shortest periodic geodesic 



(0,oo) B t 



1 (1 + it) + (-1 + u)Vs 



2 l + it 

of the modular curve PSL2(Z)\H. 

Remark 9.1. Identifying an infinitely long cylinder of circumference 1 with 
the quotient space of C under translations of the form 2i7rZ, the usual 
exponential function transforms the cylindric model into the logarithmic 
model with points on a logarithmic ontogenetic spiral defined by N 9 n 1 — > 

n n p 2iir6n 
r e 



9.2 van Iterson's disc-packing model 

In [5] van Iterson considers periodic packings of equal discs on cylinders such 
that every disc touches (at least) two pairs of adjacent discs. 

Levitov in [H] observes that the associated lattices correspond to ele- 
ments of norm 1 in the modular domain and that the corresponding Te- 
ichmuller space (given by equivalence classes of lattices endowed with a 
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positively-oriented basis) is a 3— regular tree with mid-edges given by the 
PSL2(Z)— orbit of i (corresponding to the square-lattice) and vertices given 
by the PSL^Z) — orbit of lih/l (corresponding to the hexagonal lattice). 
The rooted subtree defined in the quarter-plane of C defined by all elements 
of HI with non- negative real parts is then in natural bijection with the so- 
called Farey-tree. The optimal approximatively straight choice for a path on 
this rooted tree corresponds of course again to the golden mean and is given 
by alternating left- and right-turns at every bifurcation, as also observed by 
Levitov who gives a physical explanation based on energy levels of this fact. 

Remark 9.2. van Iterson's model involves sphere packings which are "lo- 
cally optimal" in the sense that every disc of the packing (almost) touches 
four other discs. The associated Voronoi domains have however not asymp- 
totically equal arece. 

The model determined by the phyllotactic map ip$ gives Voronoi domains 
with asymyptotically equal areasbut leads to discs in the corresponding sphere 
packing which are almost all isolated. Exceptions are occuring at parastichy 
transitions of type I (happening asymptotically at the square lattice if 6 is 
the golden ratio). 

10 Testing the existence of phyllotactic geodesies 
in real sunflowers 

Phyllotactic geodesic are perhaps a mere mathematical artefact due to the 
use of the model maps (pg. This Section sketches a test probing the reality 
of the theory. 

A first step is of course gathering real data, consisting of a fair num- 
ber of pictures of large flawless sunflower-capitula. These pictures should 
be enriched by adding as smoothly as possible (using perhaps splines or 
trigonometric functions and a least square method) all visible parastichy 
spirals. Intersections of transversal parastichy spirals should now be taken 
as the centers of seeds. Points near the center can be neglected. 

We can check adequacy of (fie(N) for sunflowers as follows: Determine 
for each picture (endowed with a complex coordinate system) parameters 
A G C, C € C*,#,7 € R giving the best least square approximation of the 
obtained seed-centers with a suitable set of points of the form 

N B n .— > A + Cy/n~+^ e 2Mn , 

supposing that the pictures have no distorsions (additional parameters are 
necessary otherwise). If this approximation is nearly perfect, the sunflower 
map is an accurate description of reality and the existence of phyllotactic 
geodesies is confirmed. A failure or a bad match does however not contradict 
the existence of phyllotactic geodesies but forces us to compute points of 
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"hypothetical" geodesies using the real data-sets instead of the model set 
¥*(N). 

This can be achieved as follows: For each point P neither on the bound- 
ary nor in the center of the sunflower, we determine pairs of points a, A 
and b, B adjacent to P with a, A on one parastichy spiral through P and 
b, B on the other, transversal parastichy spiral through P. The linearized 
lattice at P is then approximatively given by Z^y^ + Z-^^. This allows 
the computation of the corresponding modular invariant by considering the 
point of the modular curve represented by i^Ef (f° r the unique sign choice 
leading to a strictly positive imaginary part). Suitable lifts of these points 
to H should now lie close to a hyperbolic geodesic which can be guessed by 
least square approximation. 
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